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SUMMARY 

The effect of varying the design of supersonic biplanes has been 
theoretically investigated to detemiine the configuration required 
for optimum aerodynamic performance. The investigation was chiefly 
concerned with biplanes having lower and upper airfoils of equal chord 
length and of trieingular cross section. For such biplanes the changes 
in aerodynamic performance, resulting from varying the edge angles 
simultaneously and in pairs were calculated. Lift and drag coeff- 
cients were also calculated for a biplane having convex sections and 
for a configuration employing a small shock-reflecting surface in 
place of the lower airfoil. 

The theoretical aerodynamic coefficients of the biplanes inves- 
tigated are compared with those of an airfoil with diamond profile 
and with those of a thin flat plate. The variation of the center of 
pressure with angle of attack and the relative loading of the airfoils 
was also investigated for several biplanes. For one biplane of 
triangular cross section, the variation of the aerodynamic coefficients 
with flight Mach number was calculated. A discussion of the effects 
of friction drag on the relative performance of biplanes and single 
airfoils is included. 

The calculations show that, in a frictlonless supersonic air 
stream, biplanes of triangular cross section yield hi^er lift-drag 
ratios than diamond airfoils of the same thickness ratio and that, 
for high lift coefficients, unsymmetrical biplanes yield higher lift- 
drag ratios than symmetrical biplanes. When friction drag is con- 
sidered, the calculations show that biplanes with the lower airfoils 
thicker than the upper airfoils should have higher lift-drag ratios 
than symmetrical biplanes. 

For each of the biplanes an optimum spacing was found at each 
Mach number. Although the performance of the biplanes for this optimum 
spacing was found to be improved over that of a diamond airfoil, the 
calculations showed that with a constant biplane spacing this improve- 
ment was maintained over only a limited range of Mach numbers near the 
optimum. 
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INTRODUCTION 

At supersonic speeds in a frictionless fluid the aerodynamic 
shape with the least drag and the greatest lift-drag ratio is the 
thin flat plate. Practical airfoils, however, must have finite 
thickness etnd consequently will have lower lift-drag ratios than the 
thin flat plate. For a given finite thickness ratio the lowest drag 
and the greatest lift-drag ratio are obtained with airfoil sections 
having a diamond profile (fig. 1(a)). The aerodynamic performance 
of a diamond airfoil, howevei-*, may theoretically be exceeded by means 
of a proper superposition of two aiirfoils of triangulai’ cross section 
(fig. 1(b)). The possibility of using such arrangements to approx- 
imate the aerodynamic characteristics of a thj.n flat plate was sug- 
gested by Busemann In reference 1. An analysis of such biplanes was 
undertaken by Walchner in x’eference 2. As a first approximation, 
the pressures on the inner surfaces of certain properly shaped 
biplanes were shown to be mutually canceled and only the outer two 
surfaces were shown to contribute to the wave drag. These two sur- 
faces are equivalent to those of a thin flat plate . Busemann’ s 
approximation, however, aaaimios that the expansion around the inner 
corners (fig. 1(h)) taices place across a single discontinuity plane; 
whereas the expansions actually occur through a wedge-shaped region. 

A part of the expansion wave from each of the inner comers la thus 
intercepted by the rear surface of the opposite airfoil and the rest 
passes outside the biplane. (See fig. 2.) The pressures are there- 
fore not equalized internally and transition to free-stream condi- 
tions must take place externally by means of ccaaprees ion shocks. 

These transitions represent energy losses that appear as increases 
in the drag of the biplane arrangements. 

Walchner (refex'ence 2) showed that the expansion waves, can be 
completely contained within a biplane if the trail ing-odge contours 
are so curved that the rarefaction waves are not' ref lected from the 
surfaces. The required contour, which heia zero traillng-edge angles, 
iS incompatible with the strength requirement of practical wings. 
Walchner integrated the theoretical pressures over the surfaces of 
two possible biplane arrangements to determine bhe extent to which_ 
the drag and lift coefficients are altered when nonzero edge angles 
are maintained. He Concluded that when friction drag was considered 
these biplanes were approximately equal in drag to a biconvex airfoil 
of the same thickness as one of the airfoils of the biplane. 

An experimental investigation of a symmetrical supersonic biplane 
was reported by Ferri (reference 3) . The optical observations reported 
Indicate that the starting characteristics of supersonic biplanes are 
in many ways similar to those of a convergent supersonic diffuser. 

When the speed of the airstream past the biplane was increased from 
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subsonic to supersonic values, vith the spacing between the two air- 
foils set at the theoretical optimun value, the theoretical shock 
configuration was not obtained. Instead, fluctuating shock configur- 
ations were observed until the design velocity was reached; then a 
shock wave curved inward from the leading edges was observed to span 
the entrance between the two airfoils. When the spacing between the 
airfoils was sli^tly increased the expected Intersecting oblique 
shock pattern appeared and remained when the spac ing was again reduced 
to optimum. Ferri found that when the curved-shock configuration was 
present the di’ag of the biplane was about six times as great as when 
the theoretically expected shock configuration was attained, A 
probable solution to the problem of reaching optimum operating condi- 
tions is a starting rocket to accelerate the aircraft to or beyond 
its design Mach number. For wind-tunnel investigations, of course, 
the spacing may be made variable and adjusted to its optimum value 
only after the design Mach number is attained. 

The present theoretical Investigation was undertaken to deter- 
mine the effect of design variations and operating variables on the 
aerodynamic performance of biplanes. The design variations for 
triangular-section biplanes Included (l) variation of all edge angles 
simultaneously, (2) variation of the edge angles of the lower and 
upper airfoils separately, and (3) separate variation of the trailing- 
edge angles. Lift and drag coefficients were also calculated for a 
biplane having circular convex sections and for a configuration 
having a small shock-reflecting surface in place of the lower airfoil. 
The calculated aerodynamic coefficients are ccanpared with those of a 
diamond airfoil and vith a thin flat plate. The variation of the 
center of pressure with angle of attack was investigated for a 
symmetrical and an imsymmetrical biplane, and the loading cf these 
biplanes was ccaapared with that pf a diamond airfoil. For one 
biplane, the variation of aerodynamic coefficient with Mach number 
was also determined. The effect of friction drag on the relative 
performance of the biplanes and single airfoils is discussed. 


DESCRIPTION OF BIPLANS ARRANGEMENTS 

The biplanes investigated are chiefly of the type shown in fig- 
ure 1(b), that is, biplanes consisting of two triangular -section 
airfoils of equal chord length. The terras used in discussing such 
biplanes are defined in the figure. The ccrresponding terms for 
comparison with the diamond airfoil are defined in figure 1(a). Frcm 
the definition of thickness ratio, a diamond airfoil with the same 
edge angle as the thinner of the airfoils of the biplane has the same 
thickness ratio as that biplane. 
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For optimum performance, the spacing between the two airfoils 
of a triangular-section biplane (d/c)^^ is determined by the 
requirement that the shocks frcaa the leading edges shall intersect 
the surface of the opposite airfoil at the Inner turning corner. 

This optimum spacing varies with free-stream Mach number Mq and 
with leading-edge angle as shown in figure 3. The variation 

of optimum spacing with angle of attack was found to be slight and 
is not indicated on this figure. Optimum spacing was assumed through- 
out the calculations with the exception of those reported in the dis- 
cussion of the effect of varying the Mach number with constant spacing. 

Whether the biplane was rotated as a whole about a fixed axis or 
each airfoil was separately rotated about its leading edge, as in fig- 
ure 1(b), was found to be loHnaterial in determining the effect of angle 
of attack on biplane performance. The configuration with leading edges 
on a common vertical line was assumed Jn the calculations. 


METHOD OF CALCULATION 

The following symbols are used in the discussion and the figures: 
c chord ... .. 

Gj) drag coefficient for nonviscous flow 

Op ^ friction-drag coefficient 

Cj^ lift coefficient 

Cj^ pitching-moment coefficient en’ound leading edge 

Cp center-of -pressure coefficienb, e/o 

d distance between airfoils of a biplane 

d/c biplane spacing 

e distance frcsa pitching -moment axis to center of pressure 


1 



M Mach number 

P total pressure 
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p static pressiira 

t wa.YjTrmTn thickness (of thinner airfoil for "biplanes) 

t/c thickness ratio 

a angle of attack 

P Mach angle 

7 ratio of specific heats 

0 edge angle 

angle "betveen local flow direction and free-atream flow 
direction 

cp angle between shock and flow direction ahead of shock 

'P angle throu^ which flow is turned (Prandtl -Meyer theory) 

■' Subscripts : 

U upper 

L lower 

0 free-stream. 

1 leading 

2 trailing 

opt optimum values 

An analysis of the flow throu^ triangular-seCtion supersonic 
biplanes is presented in the appendix. Such an analysis shows that 
as the air passes between the airfoils it is first abruptly com- 
pressed by the deflection due to the leading edges and then expanded 
around the inner turning corners. The aerodynamic coefficients of 
the biplane are obtained by determining the pressure distribution on 
the surfaces resulting fran compression and exp^sion. The expansion 
process is readily followed with the help of the Prandtl-Meyer theory 
of flow around corners (reference 4) . This theory gives the ratio of 
static to total pressure p/P^ Mach number M^ and Mach angle p 
as functions of the angle through which the flow is turned These 

relations^ where 'P is taken eq.ual to 0 for M = 1.0^ are plotted 
In figure 4. 
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The ccmpreseive turning at the leading edges takes place throu^ 
ohliiiue shocks; the direction of these shocks is determined by the 
Mach number of the free stream and by the angle of inclination of the 
biplane STArfaoes to the free-stream direction. The pressure resulting 
from such compressive turning may be deteivained from obllq.ue-shook 
r o lat ions , The Prandt 1-Meyer ' re lat ions , however , give suf f ic ient ly 
good approximations to pressure resulting from compressive turning if 
the shocks are not too Intense; that is, if the process is approxi- 
mately isentropic. For tlie angles of inclination and the Mach numbers 
considered in the present investigation, errors resulting from the use 
of the Prandtl-Meyer theory in place of exact oblique-shock relations 
to determine pressures resulting from compressive turning were found 
to be insignificant. This tlieory was therefore utilized for conven- 
ience in calculation to determine all pressures as well as Mach lines. 
The angles between the shocks and the flow directions, Tdiich determine 
the optimum distance between the airfoils, differed considerably from 
the Mach angles assumed by the Prandtl-Meyer theory and were therefore 
obtained from exact oblique-shock relatione. 

The pressure distribution on the biplane surfaces may be obtained 
either by grapJiical or by analytical integi’atlon. Because the inte- 
gration process is rather laborious, an analytical method suitable for 
solution with computing machines was developed. This method is des- 
cribed in the appendix. 

An example of the graphical determination of pressure distribu- 
tions is given in figure 2. The angle $ Jhat determines the pres- 
sure ratio p/p, the Mach number M, and the Mach angle p is 
indicated in each region together with the_ angle between the local 
flow and the free-stream flow 7s. The continuous -expans Ion regions 
have been replaced by successions of Mach lines, eeich of which turns 
the flow through an angle of 1°. The expressions giving the inter- 
cepts of these Mach lines on the biplane surfaces in terms of local 
Mach angles are given in the appendix. 

In the calculation of the aerodynamic coefficients, the following 
simplifying assumptions were made: (l) The continuous-expeinsion 

regions can be replaced by a succession of- Mach waves, each of which 
expands the flow throu^ an angle of 1^; (2) these Mach waves are 
abruptly deflected at a definite point in the interaction region 
(see appendix); and (3) additional drags due to viscosity will be 
additive. With regard to the first two assumptions, it was felt that 
the small additional accuracy to be expected from eissumlng smaller 
expansion Intervals or plotting the Interaction region more carefully 
did not warrant the additional labor required. The third assumption 
implies that the calculated values of 'Sre correct for viscous. 
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as veil as nonviscous, air flov suid that the friction-drag coefficient 
Cj) f may t)e experimentally determined from, the difference between the 
total drag coefficient and the drag coefficient calculated for non- 
viscous flow Cjj. 

At hi^ Mach numbers the ratio of static to total pressure p/P 
becomes very small (fig. 4) and the differences between p/P on the 
front and rear or the top and hottcm, surfaces of an airfoilj which 
determine the aerodynamic coefficients, are of the order 10 for 
{4ach numbers greater than about 3.5. For high Mach numbers, there- 
fore, the accuracy of calculated aerodynamic coefficients is chiefly 
limited by the low values of p/P. For this reason, few calculations 
were mads for Mq>3.0. 

If the thickness ratio of the airfoils is very small (values 
of 0<4°), the variation of Cj, and of a supersonic biplane 

with Mach number and angle of attaci may be determined approximately 
by means of simple eq[uationa given in reference 5. These equations 
were derived on the assumptione that the expansion around the inner 
comers takes place across a single discontinuity plane, that the 
vai'letion of pressure with flow angle is linear, and that there is 
no del' lection of the compression and expansion waves at their inter- 
section points. If the biplane spacing is optimum, the equations 
reduce to those obtained for a thin flat plate. 


EFFFDT OF VARYING EDGE ANGLES 

The effect on drag and lift-drag ratio OjJC-q of varying the 
edge angles 6 of triangular-section biplanes is shown in figure 5. 
The values shown eire for the optimum spacings (d/c)Qp.j. presented 
in figure 3. Figure 5(a) shows the effect of simultaneously varying 
all edge angles. Calculations were made for three free-stream Mach 
numbers Mq and for angles of attack a of 0° and 3°. For Mq = 1.6 
the calculations were carried out only for 9 <7° because the flow 
between the airfoils hecomes subsonic \dien 0>7.3° and the 
supersonic -Dip lane theory no longer applies. The lift-drag ratio 
Cl/Cd is seen to be almost independent of Mach number for the optimum 
spacings assumed. 

The effect of varying only ©g while. 6-j^ is held constant at 
10° is shown in figure 5(b). Curves are shown for values of Mq 
of 2.0 and 3.0 at an a of 0*^ and for an Mq of 3.0 at an a of 3° 
For these conditions reaches an optimum for a of about 

5°. Lift is obtained at an a of 0° when ©g is somewhat smaller 
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than ©L "because an unsyrametrlcal flow exists between the airfoils, 
which results in a greater average pressure on the upper airfoil than 
on the lower airfoil. As ©y reduced below the optimiun, the drag 
begins to increase more rapidly than the lift and Cl/Cd again drops. 

The effect of increasing ©l while holding ©y constant is 
shown in figure 5(c) for a ©y of 7°. Curves are dra’/m for an My 
of 3.0 and for an a of 0° and For a = 3'^, Ct/Cd remains 

almost constant for ■ ©l <10°, and decreases for larger ©y. For 
a = 0°, CiyCy rapidly increases as the value of © increases to 
8.5° and then drops slightly for larger angles. 

The effect of varying the value of Sg while ©]_ la held con- 
stant at 10° is shown in figure 5(d). The optimum value in this case 
again is the result of opposing tendencies; As ©g decreases, the 
expansion around the inner comers becomes less and the average pres- 
sures on the rear inner surfaces tend to increase. The expansion 
waves from the inner corners, however, are intercepted by increasingly 
larger portions of the inner rear surfaces, and lower pressures toward 
the trailing edges result. Because the" first tendency lowers drag and 
the second tendency increases drag, an optimum value of ©2 results. 


COMPAEISOK OF P0LA3 DIAGRAMS 

For comparison with the diamond airfoil. and with the thin flat 
plate, two biplanes were chosen from figure 5 and their drag and lift 
coefficients were calculated for several additional angles of attack. 
The results are plotted in figure 6. Comparisons are made at free- 
stream Mach numbers My • of 2.0, 3.0, and 4.0 in figures 6(a), 6(b), 
and 6(c), respectively. The biplanes selected wore a symmetrical 
biplane with all edge angles © = 7° and arT unsymmetrical biplane 
with ©y = 7° ■ and ©y = 10°. The curves show that both biplanes 
give greater lift for a given drag than a diamond airfoil of the same 
thickness ratio for all values of a and My considered. For the 
biplanes, the symmetrical configuration gives greater lift for a given 
drag over.the lower range of Cy; whereas for high values of Cy the 
unsymmetrioal configuration has lower drag than all others including 
the thin flat plate. (For an My of 4.0 it was considered unnecessary 
to calculate the polar diagram for the symmetrical biplane, inasmuch as 
the relative position of the four configurations considered seems to be 
insensitive to Mach number.) 
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The ccanpai'isona presented in figure 6 indicate that the unsym- 
ffietrioal biplane has lower drag than the syuniietrlcal biplane for higjh 
values of The maximum, value of obtainable (slope of a 

strai^t line from the origin tangent to the polar curve), howevex*, 
is scmewhat less for the unsyamietrical biplane. A polar curve for a 
biplane with a ^ of 7° and a of 13° was calculated to deter- 

mine whether this maximum value of Oj^/C-q is reduced or increased 
by further increases in the values of &j_, (fig. 7). The trend toward 
greater values of Cl/Cp at hi^ values of C]^ is continued but the 
niaximum. value of Cj^/Cp is reduced as the value of 0^ is increased. 

The same three biplanes are c<aax)ared in figure 8 with diamond 
airfoils of three thickness ratios and with a thin flat plate. In 
this figui-e Cj^/Cp (which is substantially independent cf Mq) is 
plotted against a. The symmetrical biplane and the biplane with a 

of 10^ have maximum values of Cp/Cp greater than the symmetrical- 
diamond airfoil with values of 8 of 5°; that is, for a given value 
of t/c the biplanes in frictionless flow have considerably greater 
lift-drag ratios than the symmetrical-diamond airfoil. 


VAEIATION (F CENTER OF ERESSUEE 

In order to determine the center of pressure of the biplanes 
compared in figure 6, their pitching-moment coefficients about the 
midpoint between the two leading edges were calculated. The center- 
of -pressure coefficient was then obtained from the relation 
Cp = e/c = Cm/CL, where c is the chord length and e Is the dis- 
tance from the pitching -moment axis X to the center of pressure. 

(See fig. 9.) This ratio is plotted against a for the two biplanes 
of figure 6 and for the diamond airfoil of the same thickness ratio. 

The values of Cp were found to be substantially independent of Mach 
number for the biplanes as well as for the diamond airfoil. The value 
of Cp varies widely with a for the unsymmetrical biplane but is 
almost constant for the symmetrical biplane and for the diamond airfoil. 


EFFECT OF VAEYIN& MACE NUMBER WIIH CONSTANT BIPLANE SBADim 

The results presented thus far have dealt with biplanes of optimum 
spacing, which varies onlj' slightly with angle of attack but Quite 
widely with free-etream Mach number. Ferri (reference 3) found that, 
when the biplane spacing was less than or greater than optimum, the 
observed flow patterns differed greatly frcm those theoretically pre- 
dicted. The high pressure beyond the intersection of the oblique 
shocks (fig. 2) was appai’ently transmitted throiigh the boundary layer 
and resulted in a flow separation either ahead of or after the inner 
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turning corners, depending on -wiietlier the spacing vas less than or 
greater than optinium. The ohligue shocks did not strike the wing sur- 
faces hut were reflected as expansion waves from the separated fluid 
layer neai’ the surfaces. The experimentally determined lift and drag 
coefficients wei^e nevertheless found to vary with biplane spacing much 
in the manner predicted by theory (see fig. 55, reference 3) although 
the experimental variations were less than the theoretical. 

Calculations were made to determine, at least qualitatively, the 
effect of varj’-ing the value of Mq while a constant value of d/c 
was maintained. The biplane selected was the unsymmetrical one with 
©U = 7° and 6 -^ = 10°. (The analytical integration method described 
in the appendix applies only when the shocks intersect the biplane 
surfaces at the inner turning comers; consequently, the pressure dis- 
tributions for nonorrbimum spacings had to be determined graphically.) 
The variation of drag coefficient Cjj, center-of -pressure coefficient 
Op, and lift-drag ratio C]^/Cj) with f'lach nimber for this biplane is 
shown in figure 10. The biplane spacing d/c was held constant at 
a value of 0.15, which is optimum for this biplane at a value of 
Mq of 3.0. (See fig. ,3.) Calculations wero made for values of a 
of 3° and 5°. The curves show that, for these values of a, Cj^/Cj) 
remains greater than that of the diamond airfoil of the samo thickness 
ratio for a range of Mach numbers between 2.7 and 3,4. Inasmuch as 
experimental variations were less than theoretical (reference 3), the 
actual range of Mach numbers for which 0^/0^ remains greater for the 
biplane than for the diamond airfoil will probably be wider than that 
Indicated in figure 10. ■ 


EELATrVE LOADING OF BIPLANES AND DIAMOND-PEOPILE AIRFOILS 

Because the internal pressure in a biplane la considerably 
greater than atmospheric, the upper airfoil cf a biplane may be 
expected, for a given lift, to be more heavily loaded than a s Ingle 
airfoil. The relative loading is shown in figure 11. The leading 
factor A(p/P) plotted in this figure is the difference between the 
average ratio of static to total presstire on the upper and lewer 
surfaces of the airfoils. This factor is plotted eigainst lift coef- 
ficient Cj^ for a symmetrical biplane, an unsymmetrical biplane, 
and a diamond airfoil for free-stroam Mach numbers Mq of 2.0 and 
3.0. At the lower value of Mq, the tuasymanetrical biplane is more 
heavily loaded; whereas at an Mq of 3.0 the two biplanes are about 
equally loaded. Both biplanes ore more heavily loaded than the diamond 
airfoil of equal thickness ratio except at very hi^ lift coeff ioionts . 


I 
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The heavier loading of the biplanes, however, is probably no serious 
disadvantage; it should be feasible in practice to strengthen the 
biplanes by fastening the two aisrfoils together in some manner. 


EFZEXJT OF FRICTION 

Inclusion of friction effects should somewhat reduce the rela- 
tively greater lift-drag ratio of biplanes as compared with single 
airfoils, inasmuch as the additional drag due .to friction will be 
about twice as great for the biplanes. No adeq.uate data on friction- 
drag coefficients Cj) f at suporsonic speeds are yet available. The 
magnitude of the Cj) f that will reduce the lift-drag ratio of the 
biplane to that of the diamond-profile airfoil may be estimated. If 
Cd f is assumed to be twice as groat for biplanes as for single 
airfoils, the diamond airfoil will have a lift-drag ratio equal to 
that of a biplane when 


%b _ *^L,d 

%,b + *^,d + ^D,f 

or when 

„ ^L,b ‘^.d ■* ^L,d '^D.b 

- 201,^4 - 

where subscripts b and d indicate the biplane and the diamond 
airfoil, respectively. A few limiting values of Cj)^f wore calcu- 
lated from the curves of figures 6 and 7. Lift coefficient Cx, attd 
drag coefficient for frictionlese flow Cj) were calculated for the 
points of maximum Cl/Cj) for values of Mq ■ of 2.0 and 3.0, The 
results appear in the following table: 
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(dag) 


Cd 

°D,f 

Mq = 2.0 

Diamond airfoil, 
0=7° 

Biplane, 0g = 7° 
Biplane, 0 q = 7° 

7 

10 

8.15 

14.7 

14.0 

0.0164 

.0061 

.0092 

0.0037 

.0063 

Mq = 3.0 

Diamond airfoil, 

0 = 7° 1 
Biplane, 0g = 7° 
Biplane, Sy = 7° 
Biplane, 0^ = 7° 

7 

10 

13 

1 

8.15 

14.2 

13.2 
10.5 

0.0111 

.0034 i 
.0072 
. 0142 ! 

0.0017 

.0042 

.0330 


An examination of this tahlo shove that a greater value of Cj) ^ is 
allowable for low values of Mq and for largo values of These 

results indicate that, for a given value of t/o, the unsymmetry of 
the hiplano should be increased the higher the friction coefficient 
encountered . 

Similar calculations may be made to determine the limiting values 
of Cj) f above which the unsymmetrical biplane has a lift-drag ratio 
hi^er than that of the symmetrical biplane. ' Calculations using the 
same values of C^ and Op as in the table show that for an Mq 
of 2.0 the biplane with 0^ “ will yield a lift-drag ratio equal 

to that of the biplane with all edge angles 7° when . Cj> f oquads 
0.00024. At an Mq of 3.0, the biplanes with values oh 0^ ^0° 

and 13° will equal the lift-drag ratio of the symmetrical biplane when 
Cjj f is 0.00018 and 0.00041, respectively. These low values of Cjj f 
inilcato that under actual test conditions the unsymmetrical configur- 
ation will probably attain a hi^er lift-drag ratio than the symmetri- 
cal configuration. 


OTHER TiPES CF BIPLANE 

The use of circular convex sections in place of the triangular- 
section airfoils greatly increased the drag for a given thickness 
ratio t/c. For an Mq of 3.0, an a of 0°, and a t/c of 0.088, 
for example, the drag coefficient Cp with circular convex sections 
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was found to be 0.0443 as compared with. 0.0048 for the triangular- 
section biplane of the same thickness ratio. A biconvex single airfoil 
of the same thickness ratio has a Cjj of 0.0146. Hb further calcu- 
lations for cxirved-section biplanes were undertaken. 

Another method of obtaining a partial pressure cancellation con- 
sisted in replacing the lower airfoil with a small shock-reflecting 
surface, (See fig. 12.) For this arrangement, Cl and C© were 
estimated by means of graphical Integration for an Mq of 2.0, an 
a of 3°, and edge angles 0 of 7°. As expected, the results are 
dependent on the- ratio of the chord of the shock-reflecting surface 
to the chord of the airfoil L/c. The valur^s obtained are: 


L/c 

Cd 

Cp/Cp 

0.25 

0.031 

4.67 

.15 

.029 

5.58 

.10 

.028 

5.69 

.05 

.027 

5.58 


Comparison with figure 8 shows that Cl/Cjj does not approach the 
value obtaihable with a diamond -profile airfoil. The drag coeffi- 
cient, moreover, is about twice as great as for two symmetrical air- 
foils. For these reasons and because of its structural disadvantages, 
this scheme is probably of no practical interest. (For the airfoil 
of fig. 12 without the shock-reflecting surface, Cp is 0.028 and 
Ch/Cp is 5.20.) 


COWCUJSIOMS 

From calculations of the aerodynamic coefficients of biplanes 
in a frictionlese supersonic air stream, the following conclusions 
may be drawn: 

1. For any given lift coefficient, a triangular-section biplane 
has lower drag than a diamond airfoil of the same thickness ratio. 

2. When all edge angles are simultaneously varied, the lift-drag 
ratio increases as the angles are decreased, 

3. With an upper airfoil of constant thickness, the maximum lift- 
drag ratio obtainable decreases as the thickness of the lower airfoil 
is increased. For low angles of attack, however, the lift-drag ratio 
of a biplane may be increased by making the lower airfoil thicker t_ian 
the upper airfoil. If the lower-edge angles are too greatly increased, 
€in optimum value is passed and the lift-drag ratio drops. 
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4. The. lift -drag ratio of hiplanes may he slightly improved hy 

making the trailing-edge angles somewhat smaller than the leading- 
edge angles. 'r; - : - 

5, The variation of center of pressure with angle of attack is 
considerably greater for unsymmetrical than for symmetrical hiplanes. 

Inclusion of friction effects scar.ewhat modifies the foregoing 
conclusions because the additional dreg duo to friction will he about 
twice as great for biplanes as for single airfoils. Calculations 
indicate that, as friction-drag coefficient^ increases, the thickness 
ratio of the lower airfoil of a biplane should be Increased to main- 
tain lift-drag ratios higher than those of symmetrical-diamond airfoils 
and that unsymmetrical biplanes will probably yield greater lift-drag 
ratios than symmetrical biplanes in frictional flows. 

Certain practical disadvantages connected with biplanes should 
be considered in evaluating their practical usefulness. The spacing 
between the airfoils must be made variable if optimim performance is 
desired over a wide range of flight velocities. Over a certain 
limited range of flight Mach numbers, however, the biplane would 
maintain a higher lift-drag ratio than the .diamond airfoil of the . 
same thickness ratio. ..The greater loading j?f .biplanes for a given 
lift can probably be dismissed as unimportant because it should be 
feasible to increase their strength by fastening the two airfoils 
together . 


Flight Propulsion Eesearoh Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, December 10, 1946. 
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APPENDIX - DETERMINATION OB’ AERODYNAMIC 
COIEFICIEKTS OF aJPFBSQNIC BIPLANES 

Bj use of the notation of figure 13, the pressure distributions 
over the biplane surfaces may be determined as. foUovs: Upon entering 

field (l, U) the air stream is tuomed through an angle y + a) 

from the fi^ee-stream direction. The transition to field (l, U) taJces 
place througii a compression shock. The angle cp]_ n that this shock 
makes with the free-stream direction ma^ be determined from oblig.ue- 
shock relations. Similarly, upon entering field (l, L) the flow is 
turned through an angle (0^ ~ a), also in a cc^pressive ^enae, 

and the oblique shock makes an angle qp-^ with the free-stpeam 
direction. The resulting pressiires in fields (1, U) and (1, L) may 
be obtained either from oblique-shock relations or, if (6l,U + “') 
is fairly small, they may be closely approximated by the Praadtl-Moyer 
relations for flow around a comer. The Prandtl-Meyer relations are 
plotted in figure 4, Thtia, for example, if Mq is 2.0, ©■]_ y is 7°, 

and a is 3°, pj^ y/P is obtained by subtracting 10° from ihe turning 
angle 4 corresponding to an Mq of 2.0; that is, 

'^1,U = 26.3° - 10° = 16.3°; P3 y/P = 0.216. 

At the intersection of the two shocks from the leading edges, 
both shocks are deflected and the flow passes into a ccsaaon field 
(U, L) . The angle 93, U shock makes with the flow direction 

in field (l, U) is now determined from the known y and the angle 
of deflection throu^ the shock. The aiigle of deflection, in turn, 
is determined by the requirement of a coiranon flow direction in 
field (U, L) . This ozsiEiLon direction can be attained only if the sum 
of the deflections through the two upper shocks equals the sum of the 
deflections through the two lower shocks. The deflection through the 
second upper shock is therefore (0^ l "O') degi’ses and throu^ the 
second lower shock is y + a) degrees. The flow in field (U, L) 

has been turned (0q u f a + 0-]_ y - a) degrees in a compressive sense 
and the conditions in field (U, L) can now be road from figure 4. 
Aithou^ the turning is compressive through both shocks, its direction 
is reversed and the actual flow direction in field (U, L) is 
^U,L = (61, U t a) - (0i^L - tx) = (01^XJ - ®1,L + 2a) degrees from the 
free-stream direction. 

In passing into field (2, U) the flew is expanded 
(®2,U ' ^ ^JJ degrees from its direction in field (U, L) . The 

flow in field (2, l) is expanded (0g y a - degrees. The 

expansions take place through' a wedge-s^ped region with limits 
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determined by the values of p corresponding to the conditions before 
and after the expansion. The conditions are again obtainable fi-om 
figure i because the turning angles are known. 

The expansion waves are mutually deflected in the region of inter- 
section. The angles that the bounding Mach lines beyond the inter- 
action zone make with the flow direction are again determined by the 
conditions before and after the expansions. The angle of the line 

closing off field (2^ U)j for example^ is |sin”^(l/M 2 y) f (6g y - ot.)J . 
The angle of the line terminating the upper expansion region is deter- 
mined by the number of degrees from field (2, U) that the flow is 
turned by the expansion wave. Each expansion wave is reflected from 
the surface and passes outside the airfoils. The four shocks ema- 
nating frcmi the two trailing edges bring the entire flow pattern back 
to free-stream conditions at infinity. The angles of these shocks 
can be determined from this requirement (that is, free-stream condi- 
tions at infinity), but this detoi’mination is unnecessary to find the 
aerodynamic coefficients, for conditions beyond the trailing edges 
have no effect on the biplane. 

In order to integrate the pressures over the airfoil surfaces, 
the continuous expansion waves are replaced with a succession of .M^h 
waves, each of which turns the flow through a small angle &9, 
Throu^out this investigation the value A0 was fixed at 1°. The 
problem of integration is obviously simple except over the portions 
of the inner rear surfaces that intercept the expansion waves from 
the opposite airfoil. This region may be treated analytically, by 
obtaining expressions giving the point at which each of the succession 
of Mach waves intersects the opposite surface. A suimaation process 
can then be made to determine the drag and lift forces contributed 
by these regions. 


The relations between the angles of the Mach lines and their 
intercepts on the airfoil surfaces may be deduced from the sine law. 
The assumption is made that the bending of each Mach lino through the 
interaction region takes place .entirely at the point where it inter- 
sects the corresponding Mach line emanating from the opposite airfoil. 
The following expressions are obtained by using the notation of fig- 
ure 14 : 


°1,U sin (9i,u " ^l,U ~ ”') sin - a) 

^ " cos 01, u sin •(- - ®1,U " 

_ °1,L sin (cpp^x, - 01 ^ 1 ^ + g) sin (cpg^j;^ - h- a) 

^ ^1,L ('^,L ^ “ ®1,L + 


( 2 ) 
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L t^ Pj^ n 

(3) 


1 4 - 1 ■ 

tan Pi^L 



tan Pq Tj 

(4) 


Ftt - ^-r 

U -L tan pq^L 

Ly = 

p (^1,U + ^2,U) 

/ C \ 

U sin Pq^u sin [^2,U ~ ®2,U 

w; 

It. = 

sin (Pq^q + 

(6) 

L sin pq T sin (Pg q - ®2,L “ 


where 

D distance frcm shock intaraectlon to Inner turning comer 

F distance -from intersection of expulsion waves to inner turning 

comer 

L. distance from inner turning comer to intersection of expansion 
vrave with biplane surface 

In each case p is the angle that the Mach lines make with the free- 
stream direction; that is^ it is the Mach angle plus the angle which 
the local flow makes with the free-stream direction. Because the 
flow ‘direction changes 1° through each Mach line, the expansion process 
must be followed (juite ceurefuUy to determine the correct values of p 
to use in equations (o) to (6). When more Mach lines ememate from one 
of the airfoils, supplementary Mach lines may be assumed to emanate 
frcm the opposite airfoil in order that the bending points of the 
excess lines may be determined from equations (3) to (6). 

When the intercepts of the Mach lines on the airfoil surfaces are 
known, the aorodynamic coefficients per unit span are obtained by a 
summation of pressures OTer the surfaces. 

The course of the expansion may be clarified by studying the 
example of figure 2. Here the values of 'H to bo used in figure 4 
and the direction of the flow A are Indicated in each region. The 
conditions I-Iq = 2.0, 6^ tj = ®2,L = •> 

a. = 3° were assumed for this sketch. 
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The analytical method Just described is useful only when the 
shocks from the leading edges intersect the opposite airfoil at the 
inner turning comer. If the coefficients for other cases are 
desired, a graphical integration must be made or some other analytical 
expressions derived. Graphical integration was used in the present 
investigation to determine the effect of varying the Mach number while 
a fixed biplane spacing was maintained and to determine the coeffi- 
cients for the convex- section biplane and for the single airfoil with 
shock-reflecting surface. 
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®1,U “ ®1,L = ®1 

®1,U “ ®1,L “ 

®1,U = '7°* ©l,L “ ®1 


Leadlng-*edge angle, 6|, deg 

Figure 3. - Variation of optimum biplane spacing with leading 
edge angles for several tfach numbers. 
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Drag coefficient, Cq 

Figure 7,— Effect of Increasing ttie lower-edge angles on polar 
diagram of triangular-section biplanes, Pree-stream Mach 
number Mq, 3,0j upper-edge angle 9^, 7°, 
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Mach number Mq, 3.0. 




Center>of-p res sure coefficient 


NACA TN No. J3J6 


Fig. 9 
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Figure 9.- Effect of angle of attack on center of pressure of 
two biplanes and a diamond airfoil. (Op = e/c.) Free- 
stream Mach number Mq, 3«0« 




Drag coefficient, Cjj Center-of-pressure Lift-drag ratio, 

coefficient, C_ 


Fig. 10 
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Free-strean Mach number, Mq 

' Figure 10* Variation of drag ooefflolent, center of pressure, 
and lift-drag ratio with Mach number for biplane with 
optlmun spacing for free-strean Mach number M of 3,0, 
(d/c-0,15.) 0 
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Fig. II 



Figure 11.- Relation between lift coefficient and loading 
factor for upper airfoils of two biplanes and for a diamond 
airfoil at free-^stream Mach, numbers Mq of 2.0 and 3»0» 
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Plgupe 12,- Sketch of single triangular airfoil with shock 
reflection surface, Pree-stream Mach number M«. 2.0. 
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Figure 13,— Notation for analysis of flow through biplanes 
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Fig. 14 



Figure 14*'- Notation used for analytical Integration process. 




